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Spacecraft  Charging 
at  Geosynchronous  Orbit — 
Solution  for  Eclipse  Passage 


1.  INTRODUCTION 

The  interactions  between  spacecraft  and  the  ambient  space  environment  have 
become  an  increasing  problem  to  the  space  physics  community.  The  best  known  of 
these  environmental  interactions  is  that  resulting  from  spacecraft  charge  buildup 
(DeForest1).  In  this  report,  simple  plasma  probe  theory  will  be  used  to  develop 
equations  which  allow  estimations  of  the  potential  that  exists  between  a spacecraft 
and  the  ambient  medium  in  the  presence  of  a varying  photoelectron  flux.  The  vari- 
ations are  compared  with  observations  from  the  geosynchronous  satellites  ATS-5 
and  ATS-6. 

2.  BASIC  EQUATIONS 

2 

The  basic  equation  describing  the  spacecraft  charging  phenomenon  is  (WhippJe, 
and  DeForest*): 


(Received  for  publication  11  May  1978) 

1.  DeForest,  S.  E.  (1972)  Spacecraft  charging  at  synchronous  orbit,  J.  Geophys. 

Res.  77:651-659. 

2.  Whipple,  E.C.  (1965)  The  Equilibrium  Electric  Potential  of  a Body  in  the  Upper 

Atmosphere,  NASA  X-61 5-65-296. 


(1) 


' (JI  + JSE 


+ JSI  + J 


BSE 


+ JPH>  = 


where 

J_,  = Incident  electron  current 

ht 

Jj  = Incident  ion  current 

J,,j  = Secondary  electron  current  due  to  ions 

J__  = Secondary  electron  current  due  to  electrons 
or* 

J-DOE'  = Back-scattered  electron  current 
Jpjj  = Photoelectron  current. 

To  first  order,  for  the  range  of  currents  and  potentials  observed  at  geosynch- 
ronous orbit,  the  secondary  and  backscattered  currents  are  approximated  by  linear 
relations: 


JSE  ~ a JE 

JBSE  * b JE 

JSI  * C JI 
such  that: 


(2) 


x JE  -yJi  - JpH  - 0 


(3) 


for: 


x = 1 - a - b 
y = 1 + c . 

3 

From  probe  theory  (Chen  ) for  a thick  sheath  around  a spherical  probe  or  in 
the  absence  of  a sheath  (DeForest1),  the  current  to  a probe  as  a function  of  the 
ambient  current  and  spacecraft  potential  for  a negatively  charged  spacecraft  is: 


3.  Chen,  F.  F.  (1965)  Electric  Probes.  Plasma  Diagnostic  Techniques. 
Huddlestone  and  S.  L.  Leonard  (Ed. ),  Academic  Press,  Inc.,  New 
pp  113-200. 
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R.  H. 
York, 


(4) 


JE(V)  J 


EO 


-|qv|/TE 


Jj(V) 


'IO 


/ lsxi  \ 

(1+  Tr  ) 


(5) 


where 


V = Spacecraft  potential  in  volts 

Te  = Maxwellian  temperature  of  electrons  in  electron  volts 

Tj  - Maxwellian  temperature  of  ions  in  electron  volts 

q = Unit  charge  such  that  qV  is  in  electron  volts 

= Ambient  electron  current 
EO 

Jj0  = Ambient  ion  current. 


Eq.  (3)  becomes: 

-|qv|  /t 


x JEO  e 


- y J 


I qv  | 


io 


(lt  tT  )’ 


■’PH'0 


(6) 


Jp^j,  the  photoelectron  current,  can  be  expressed  as  a function  of  the  pre- 
eclipse  photoelectron  current,  JpQ,  and  an  attenuation  function  f representing  the 
fraction  of  sunlight  falling  on  the  spacecraft.  It  is  convenient  to  express  this  frac- 
tion in  terms  of  the  minimum  altitude,  Xm>  for  the  ray  between  the  satellite  and 
the  center  of  the  sun  (Figure  1).  Thus  Eq.  (6)  becomes: 


I qV  I . 

Ti 


Ln[(y/x)(JIO/JEO)(l  + IqVl/Tj)  + (l/x)Upo/JEO)  f(Xm)l  . (7) 


Equation  (7)  is  the  basic  spacecraft  charging  equation  in  the  presence  of  photo- 
electron  emission.  Assumptions  will  be  made  that  reduce  Eq.  (7)  to  an  analytic 
expression  giving  the  spacecraft  potential  as  a function  of  the  electron  temperature. 
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Figure  1.  The  Altitude  Xm  Above  the  Earth  of  the  Minimum  Ray  Path  Between 
the  Satellite  and  the  Center  of  the  Sun 


3.  POTENTIAL  IN  ECLIPSE 


When  a spacecraft  is  in  the  earth's  shadow,  the  photoelectron  flux  approaches 
zero  (scattered  light,  moonlight,  etc.,  may  still  be  present)  so  that  Eq.  (7) 
becomes: 


< 


= - Ln[(y/x)(JIO/JEO)(l  + iqVj/Tj)] 


(8) 


Te.  Tp  and  were  calculated  for  21  electron  and  ion  spectra  from  ATS-5 

and  for  four  spectra  from  ATS-G  for  intervals  immediately  preceding  eclipse  entry 
and  immediately  following  exit  from  eclipse.  These  spectra  cover  the  energy 
ranges  50  eV  to  50  keV  and  1 eV  to  80  keV  respectively  (see  DeForest  and 

4 5 

Mcllwain,  and  Mauk  and  Mcllwain  for  details). 


f 

i 


Parameters  x and  y for  the  range  of  observed  parameters  are  estimated  from 
formulas  given  in  Whipple2  and  Garrett*’  to  be: 


4.  DeForest,  S.  E. , and  Mcllwain,  C.  E.  (1971)  Plasma  clouds  in  tin  magneto- 

sphere, J.  Geophys.  Res.  76:3587-361 1. 

5.  Mauk,  B.  H. , and  Mcllwain,  C.  E.  (1975)  ATS-6  UCSD  auroral  particles 

experiment,  IEEE  Trans.  Aerospace  and  Electronic  Systems,  AES-1  l(No.  6): 
1125-1130. 

6.  Garrett,  H.B.  (1978)  Spacecraft  Potential  Calculations  — A Model,  AFGL- 

TR-78-01 16. 
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a « 0.  4 


(9) 


I 

t 

> 


b = 0.  2 

c * 3 


Therefore: 


x = 0.  4 


(10) 


y = 4 . 

For  the  eclipse  potentials  corresponding  to  the  25  intervals,  the  ratio  |qVoj/Tj  is 
in  general  less  than  1.  Therefore,  Ln(l  + ! q | /Tj)  is  approximately  0.  Equation  8 
reduces  to: 


Iqv 


T = Ln  K 1 / 10)<JEO/JIO>]  = Ln((l/10)(nE/ni)(TE 


(11) 


where 

nE  ' Ambient  electron  density 
nj  = Ambient  ion  density 
mE  - Mass  of  electron 
m.  = Mass  of  ion 


For  the  data  studied,  T,  ~ 2.  5 T„  and  n_  = n.  so  that: 

I E El 

JEO/JIO  ' (nE/nI)(TE/TI)1/2(mI/mE)1/2  * 25  • <12) 

For  the  25  events  studied,  the  ratio  *3  actually  24  ± 7.  As  the  current 

ratio  enters  as  a logarithm,  Eq.  (11)  becomes: 


!qv0l  ' Te  En(2.  5)  * Te  . (13) 

In  Figure  2,  the  electron  temperature  (assumed  to  be  Maxwellian)  is  plotted 
versus  the  observed  potential.  An  expression  linear  in  the  ambient  temperature 
was  fit  to  the  data  by  the  least  squares  method  (the  two  marked  points  were  not 
included  as  the  plasma  is  believed  to  have  changed  between  the  measurement  of 
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SATELLITE  POTENTIAL  ( kV ) IN  ECLIPSE 

Figure  2.  The  Observed  Temperature  of  the  Ambient  Electrons  vs  the 
Satellite  Potential  in  Eclipse.  TQ(AVG)  is  given  by  2/3 's  of  the  ratio  of 
the  energy  density  to  the  number  density.  T (RMS)  is  given  by  1/2  of  the 
ratio  of  the  energy  flux  to  the  number  flux.  The  difference  between  these 
two  temperatures  is  a measure  of  the  adequacy  of  a Maxwellian  approxi- 
approximation 


Tj.  and  the  potential).  Depending  on  which  definition  of  the  temperature  is  used  the 
fits  and  correlation  coefficients  are: 

V * = 1391  - 1.  63  T„(AVG)  (r  = 0.  86) 

o E 

(14) 

V * = 1828  - 1.  29  TJRMS)  <r  = 0.  92) 

O L 

These  relationships  not  only  demonstrate  the  validity  of  Eq.  (13),  but  imply  the 
existence  of  a threshold  electron  temperature  above  which  the  satellite  potential  is 
linear.  Such  behavior,  though  not  explicit  in  our  simple  theory,  is  due  to  the  fact 
that  at  an  intermediate  energy  (usually  a few  hundred  eV),  the  secondary  yield  is 
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greater  than  1 (Rubin  et  al  ).  The  electron  temperature  must  be  several  times 
greater  than  this  threshold  energy  before  significant  charge  buildup  occurs. 

4.  POTENTIAL  DURING  ECLIPSE  PASSAGE 

As  a spacecraft  enters  and  exits  the  earth's  penumbra,  it  sees  a varying 
illumination  generating  a varying  photoelectron  current  as  the  sun'3  disc  is  ob- 
scured by  the  limb  of  the  earth.  The  percentage  of  solar  illumination  as  a function 
of  X seen  by  a satellite  at  geosynchronous  orbit  is  plotted  in  Figure  3 (the  atten- 
uation function  is  based  on  the  work  of  Garrett7 8 9  and,  as  the  details  are  not  of  con- 
cern to  this  generalized  treatment,  the  reader  is  referred  to  that  paper).  Dividing 
Eq.  (7)  by  Eq.  (8): 

V(Xm>  = ln[(y/x)(J10/JE0)(l  + IqVl/Tj)  + l/x(Jpo/JEn)  f(Xm>l 

Vo  Ln((y/x)(JIO  /JEO)(1+  lqVol/Tl>l 

Estimating  as  before  and,  as  | q V | < I qV  | : 

V(X) 

-y2-  « - Ln[.4  + 2.5  (JpQ/JEO)  f(X  » (16) 

o 

The  ratio  J DrJJvr.  is  not  known  accurately  as  it  is  a function  of  the  material 

* vJ  ijU  2 8 1 

and  structure  of  the  spacecraft.  From  Whipple,  Grard  et  al,  DeForest,  1 and 

9 

Garrett,  a reasonable  range  of  values  is: 

^I-Jpo/Jeo-16  (17) 


7.  Rubin,  A.G. , Rothwell,  P.  L. , and  Yates,  G.  K.  (1978)  Reduction  of  spacecraft 

charging  using  highly  emissive  surface  materials,  Proc.  1978  Symposium  on 
the  Effect  of  the  Ionosphere  on  Space  and  Terrestrial  Systems,  Naval 
Research  Lab. , Washington  b.  C. 

8.  Grard,  R.J.L. , Knott,  K. , and  Pedersen,  A.  ( 1973)  The  influence  of  photo- 

electron and  secondary  electron  emission  on  electric  field  measurements  in 
the  magnetosphere  and  solar  wind,  Photon  and  Particle  Interactions  with 
Surfaces  in  Space,  R.J.C.  Grard  (Ed. ),  D.  Reidel  Publishing  Co. , Dordrecht, 
Holland,  pp.  163-189. 

9.  Garrett,  H.B.  (1978)  Effects  of  a Time-Varying  Photoelectron  Flux  on  Space- 

craft  Potential.  AFGL -TK -78-01 19, 
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In  Figure  3,  Eq.  ( 1 f»)  is  plotted  for  f(Xm>  given  by  Garrett  and  for  JpQ/JjrfJ  1.6 
and  lfi.  Also  shown  are  data  from  ATS-5  and  ATS-6  for  the  25  eclipse  passages. 
The  observations  are  subject  to  a large  (~25  km)  error  in  Xm,  as  the  orbits  of 
ATS-5  and  ATS-fi  were  not  updated  sufficiently  often  to  allow  a more  accurate 

9 

position  determination  (Garrett'  ).  Even  so,  the  predicted  variation  in  parameters 
brackets  the  actual  observations. 


Figure  3.  The  Ratio  of  the  Varying  Potential  to  the  Max- 
imum Potential  When  the  Spacecraft  Were  in  Eclipse. 
Also  shown  are  the  percent  illumination  (same  scale  as 
* V(Xm)/VJ  and  the  predicted  curves  for  IDO/Ieo  values  of 

I.  6 and  lfi 


5.  DISCUSSION 


In  the  previous  sections,  we  have  presented  equations  (Eq.  (13)  and  (16))  that 
approximate  ATS-5  and  ATS-6  observations  of  the  varying  potential  that  exists 
between  the  spacecraft  and  the  ambient  plasma.  We  believe  that  this  is  primarily 
the  result  of  the  rough  constancy  of  the  ratio  ^iq/^eo  tlie  lacl  tllat  variations 
in  this  ratio  enter  as  a logarithm.  Further,  it  is  tacitly  assumed  that  Tj  and  TE, 
the  Maxwellian  temperatures,  have  meaning.  Although  these  assumptions  are 
not  valid  at  all  times  at  geosynchronous  orbit  (particularly  at  quiet  times  when 
non -Maxwellian  ring  currents  are  observed),  they  are  adequate  following  plasma 
injections  when  a satellite  is  near  midnight.  As  this  is  a prerequisite  for  charging 
to  occur  during  eclipse  passage.  Eqs.  (13)  and  (16)  are  valid. 

Equations  (13)  and  (16)  are  of  use  in  many  space  physics  problems.  As  an 
example,  given  T£,  the  potential  in  eclipse  can  be  estimated  by  Eq.  (13).  Equiva- 
lently we  have  plotted  the  frequency  of  occurrence  of  potentials  in  eclipse  for 
ATS-5  and  ATS-6  for  4 years.  With  Figure  4 and  Eq.  (16),  the  occurrence 
frequency  of  various  potential  levels  as  the  spacecraft  entered  into  eclipse  can  be 
estimated  for  the  same  4 year  period  (Rubin  and  Garrett1®).  As  ATS-5  is  a spinning 
cylindrical  shape  (1.  3 m in  diameter  by  2 m in  length)  and  ATS-6  is  a large  (10  m) 
spin-stabilized  dish  antenna,  the  equations  are  also  applicable  to  a variety  of 
problems  and  simple  space  structures.  For  instance,  by  equating  the  photoelectron 
current  to  an  electron  beam,  Eq.  (13)  and  (16)  can  be  used  to  simulate  the  effects 
of  a varying  electron  beam. 

Although  the  scaling  to  large  space  structures  may  not  be  completely  valid, 
the  results  of  Eq.  (13)  and  (16)  may  be  used  in  principle  to  estimate  the  potential 
that  would  exist  between  the  ends  of  a large  space  structure.  This  can  be  most 
readily  visualized  by  assuming  the  structure  to  be  approximated  by  two  spheres 
separated  by  a distance  equal  to  that  of  the  structure,  and  connected  by  a long 
cable  having  the  electrical  loading  characteristics  of  the  structure.  For  infinite 
resistance,  Eqs.  (13)  and  (16)  can  be  used  to  estimate  the  potential  on  each  sphere 
as  it  passes  into  eclipse.  Assuming  the  distance  between  the  two  spheres  to  be 
equal  to  their  difference  in  X , the  potential  difference  between  the  two  ends  can 
be  estimated  as  they  pass  into  eclipse.  Depending  on  the  exact  nature  and  size  of 
the  structure,  this  simple  analogy  would  predict  potentials  as  high  as  -20,  000  V, 
a clear  hazard  to  the  spacecraft. 


10.  Rubin,  A.G.,  and  Garrett,  H.  B.  (1978)  A Statistical  Study  of  Eclipse  Poten- 
tials for  ATS-5  and  ATS-6.  to  appear  as  AFGL  In-House  Report. 
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In  conclusion,  elementary  probe  theory  suggests  basic  relations  between  the 
electron  temperature  (assumed  Maxwellian)  and  the  spacecraft  potential  relative 
to  the  ambient  plasma  when  the  spacecraft  is  in  the  earth's  shadow.  Relations 
describing  the  potential  variations  as  the  spacecraft  passes  into  eclipse  are  derived 
Observations  from  ATS-5  and  ATS-G  support  these  derivations  and  indicate  that  the 
reason  for  agreement,  is  because  the  ratio  of  the  ambient  ion  to  electron  current  is 
roughly  constant  for  injection  events  and  enters  as  a logarithm.  Provided  that  the 
results  are  of  general  applicability,  they  can  be  used  to  estimate  the  potentials  on 
large  space  structures  and  the  effects  of  electron  beams. 


t 

$ 

Figure  4.  The  Occurrence  Frequency  of  Eclipse  Potentials 
• for  ATS-5  and  ATS-G.  10  minute  averages  were  employed 
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